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Abstract 

We perform the two-loop matching calculation for heavy-to-light currents from 
QCD onto soft-collinear effective theory for the complete set of Dirac structures. 
The newly obtained matching coefficients enter several phenomenological applica- 
tions, of which we discuss heavy-to-light form factor ratios and exclusive radia- 
tive decays, as well as the semi-inclusive decay B — )• Xsi'^i~- For this decay, 
we observe a significant shift of the forward-backward asymmetry zero and find 
g2 = (3.34J:o:25) GeV^ for an invariant mass cut m'j^^ = 2.0 GeV. 



1 Introduction 

The flavour-changing quark currents qTib, with Fj = {l,75,7^,757^,io"'^'^}, govern the 
hadronic dynamics in semi-leptonic and radiative B decays. The matrix elements of the 
currents, usually parameterized by several transition form factors, are also important 
inputs to the factorization formulae for non-leptonic B decays [1]. In the kinematic 
region where the hadronic final state has small invariant mass but large energy, soft- 
collinear effective theory (SCET) |2l[3] is the appropriate theoretical framework, with 
which transparent factorization formulae for the heavy-to-light form factors have been 
derived |1] (see also [5|5] ). Thus, the accurate representation of the heavy-to-light currents 
in SCET is of particular interest. 

The LO and NLO matching coefficients for heavy-to-light currents from QCD onto 
SCET for an arbitrary Dirac matrix has been worked out a few years ago [SlIZlll]. The 
coefficients for V-A currents have recently been determined to NNLO in the context 
of inclusive semi-leptonic B decays [SHI2] in the shape-function region. In this paper 
we complete the NNLO calculation by computing the remaining matching coefficients 
of the tensor currents. The tensor matching coefficients enter several phenomenological 
applications, of which we shall discuss heavy-to-light form factor ratios and exclusive 
radiative decays, as well as the semi-inclusive decay B — )■ Xgi^i^. 

The paper is organized as follows. In Section |2] we first set up notation and then 
briefly recapitulate the techniques applied and the necessary ingredients for the two- 
loop calculation. In Section [3] the two- loop calculation of the QCD form factors and 
the corresponding matching coefficients are presented in detail. Three interesting phe- 
nomenological applications of our results to heavy-to-light form factor ratios, exclusive 
radiative decays, as well as the inclusive decay B — )• Xgi^i' are discussed in Sections S] 
and O We conclude in Section [61 The lengthy analytic expressions for the coefficient 
functions can be found in Appendices A and B. 

2 NNLO calculation 

2.1 Set-up of the matching calculation 

A generic heavy-to-light current gFj 6 is represented in SCET by a set of non-local "two- 
body" and "three-body" [SHB] operators, 

[qT,b]{0) = Y, [dsCi{s)[^Whc]{sn+)^hM 

where h^ is the static heavy quark field defined in HQET, whereas S, and Whc are the 
hard-collinear light quark field and a hard-collinear Wilson line from SCET, respectively. 
In this paper we are concerned with the calculation of the matching coefficients in the 



first line of ([T]). Tlie tliree-body operators 01 {si,S2) in tlie second line are l/irih- 
suppressed but relevant at leading power for exclusive transitions and form factors due 
to the matrix element suppression of the leading term. Their one-loop matching coef- 
ficients are known from [8l|13] and this suffices to work out their contribution to the 
exclusive transitions at (9(a^). We refer to [H] for the details of the calculation of these 
spectator-scattering terms. In the current work we consider the missing C?(af) matching 
coefficients of the two-body operators [^M^hc] {sn^)T'jhy{0) and adopt the momentum 
space representation, which follows from 



CRn+p) 



ds e'^^'+'P CUs) 



(2) 



We decompose the heavy-to-light currents in the basis from [H] (summarized in Table [T]) 
with two reference vectors v and n_ that fulfill v = (n_ + n^)/2, n^ = and n+n^ = 2. 
The matching calculation involves 12 coefficient functions Cj, which are not independent 
in a renormalization scheme with anti-commuting 75 due to the chiral symmetry of QCD. 
In the NDR scheme adopted in this work, this translates into the constraints Cp = Cs 
and C\ = Cy, while similar relations hold between the matching coefficients of the tensor 
and the pseudotensor current. As the latter is reducible in four space-time dimensions, 
we obtain the additional constraints C^ = C^ = in four dimensions. We nevertheless 
keep the more general basis from Table [T], since we work in dimensional regularization 
and obtain intermediate results that are valid in c? = 4 — 2e dimensions, where C^ and 
C^ are of 0{e) but non- vanishing. 

It is convenient to perform the matching calculation with on-shell quarks and to 
use dimensional regularization to regularize ultraviolet (UV) and infrared (IR) singular- 
ities. The SCET diagrams are then scaleless and vanish and the computation essentially 
amounts to a two-loop calculation in QCD. We, in particular, introduce an analogous 
tensor decomposition to ([1]) and parameterize the QCD result in terms of 12 form fac- 
tors, 



{q{p)mb\b{pb)) = J2 ^i^) u{v)'^',u{vb) 



(3) 



where ph = f^^bV is the momentum of the heavy quark, p = uin[,n^/2 the momentum of 
the light quark and q^ = {pb — pY = (1 ~ u)m1 denotes the momentum transfer. Due to 
the absence of loop contributions on the effective theory side, the SCET matrix elements 
are given by the tree level matrix elements multiplied by a universal renormalization 
factor Zj of the SCET current [^VF/ic] ^'jhv There is thus a one-to-one correspondence 



r. 


1 


75 


^M 


757^" 


icr''^ 


r;- 


1 


75 


r 


v^ 


nt 


757^ 


v^l^ 


^-75 


^[m^'^I 


yil^Y] 


n^^Y^ 


nK^^ 


ci 


Cs 


Cp 




fi2 


CI 


c\ 


c\ 


c\ 




Ot^ 


\_Jn-> 


C^ 



Table 1: Matching coefficients Cj according to the decomposition in ^^ (a^'^b'^' = a^b'^—a'^b^). 
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between the matching coefficients Cl and the form factors F/ 

CI = Z-'Fi. (4) 

As the form factors are, however, in general IR-divergent, there exists no analogous 
relation on the form factor level to the four- dimensional constraints C|. = C^ = 0. 

The purpose of our analysis consists in the computation of the matching coefficients 
Cj (and the respective form factors F/) to NNLO in QCD. Whereas the NLO corrections 
have been worked out in [2|l7f8] . the coefficients Cy and C\ have recently been determined 
to NNLO in the context of inclusive semi-leptonic B decays |9l4T2]. In the current work 
we complete the NNLO calculation by computing the remaining matching coefficients Cs, 
Cp and C^. The four-dimensional constraints mentioned above, will serve as a non-trivial 
check of our calculation. 

2.2 Technical aspects of the calculation 

We organize the calculation along the strategy that we used in our previous works on 
the V-A current [m[T2] . The calculation is based on an automated reduction algorithm, 
which uses integration-by-parts techniques [H] and the Laporta algorithm [TB] to express 
the two- loop diagrams (shown in Figure 1 of [TT]) in terms of a small set of scalar master 
integrals. The required master integrals are already known from the computations in [9]- 

Our results will be given in terms of the following set of harmonic polylogarithms 
(HPLs) tl9j, 

i^(0,0,l;x) = Li3(x), 
/7(0,l,l;a;) = Si,2(a;), 
H{0,0,0,l;x) = Li4(x), 
i7(0,0,l,l;x) = S2,2(x), 
i7(0,l,l,l;x) = Si,3(x), 
H{-l,0,l;x) = niix), H{0,-l,0,l;x) = n2ix), (5) 

where we introduced a shorthand notation for the last two HPLs. Whereas the ffist one 
can be written in a compact form [20] . 
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7^1 (x) = ln(l + a;)Li2(a;) + ^81,2(0;^) - 81,2(0;) - Si,2(-a;), (6) 

the second one, "^2(2;) = Jq dx'T-Li{x') / x' , cannot be expressed in terms of Nielsen Poly- 
logarithms and has to be evaluated numerically. 

The charm quark enters the matching calculation at the two-loop level through the 
gluon self energy which contains closed fermion loops. Our analytical results from Sec- 
tions 13.11 and 13.21 are valid for massless charm quark, but we also show some numerical 



results in Section 13.21 that include charm mass effects. In this case we formally keep 
mc/mf, fixed in the heavy-quark expansion, so the coefficients depend non-trivially on 
the quark mass ratio (see Section 5 of |.12j). 

The pure two-loop calculation yields bare form factors F/ that are UV- and IR- 
divergent. The UV-divergences are subtracted in a standard renormalization procedure, 
which has been described in detail in our previous works pTl[T2] . We, in particular, renor- 
malize the strong coupling constant in the MS-scheme, whereas the quark wave-functions 
and the 5-quark mass are renormalized in the on-shell scheme. The only difference in the 
current calculation consists in the fact that the scalar and the tensor current have non- 
vanishing anomalous dimensions in contrast to the vector current considered in pT|[T^. 
This gives rise to an additional multiplicative counterterm Z^^ {i = S,T). We expand 
the inverse 




Kfc) 
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in terms of the renormalized coupling constant of a theory with five active quark flavours. 
In the MS-scheme the respective NLO coefficients are then given by Z^ = SCp/e and 
Zrp = —Cp/e for the scalar and the tensor current, respectively. At NNLO the counter- 
terms can be inferred from 
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where ri/ = 5 denotes the number of active quark flavours. 



3 Results 

3.1 Renormalized form factors 

We first present our results for the renormalized form factors F/ , which are U V-finite but 
IR-divergent. It will be convenient to decompose the form factors according to 
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In this normalization the form factors become at tree level 
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(10) 



Here and below we do not quote our results for the pseudoscalar and the (axial) vector 
current, since the former are equal to those of the scalar current in the NDR scheme, 
while the latter have already been computed before and can be found in |9HT2]. 

One-loop form factors. At NLO we compute the form factors up to terms of O(e^). 
Our results are given in terms of a set of coefficient functions gi{u), which we list in 
Appendix El The scalar form factor is IR-divergent and becomes (with g^ = uml, 
u = 1 — u and L = In fi^ /ml), 
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The first tensor form factor is also IR-divergent and given by 
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(12) 



whereas the other tensor form factors are IR-finite at NLO and read 
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Two-loop form factors. At NNLO the IR-divergent parts of the form factors can be 
expressed in terms of the one- loop coefficient functions gi{u). The divergent terms of the 
scalar form factor read 



F^%{u) = C'f(l- goiu)\ + ^CaCf - uiTfCf, 
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and for the first tensor form factor we get 
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The IR-divergent parts of the other tensor form factors are given by 

F^fHiu) = -C^griu), 
and 

Ft,-i{'>J') = CFl9o(,u)gg{u) - gio{u) - 2gg{u)L 
and 
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The finite parts of the two-loop form factors involve a new set of coefficient functions 
hi{u), which we specify in Appendix iBl We find 
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3.2 Matching coefficients 

The matching coefficients C/ follow from the above expressions for the renormalized form 
factors F/ after multiplication with the inverse of the renormalization factor of the SCET 
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current Zj, cf. (jl]). To this end one has to keep in mind that the form factors have been 
computed in QCD with five active quark flavours, while Zj is usually given in SCET 
with four active flavours. We thus have 
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where n^ = n^ — 1 = 4 is the number of active quark flavours in the effective theory. 

We now expand the matching coefficients in terms of the coupling constant of the 
four-flavour theory as 
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and rewrite (j4]) up to NNLO, which yields 
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Notice that the last relation implies a term which stems from the conversion of the flve- 
flavour to the four-flavour coupling constant, 
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with (see also [II1[I2] for further details) 
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At LO the matching coefficients then become 
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At NLO the matching coefficients are given by the finite terms of the one-loop form 
factors, 
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four- dimensional constraints for the tensor coefficients that we mentioned in Section 12.11 
Here and in the following we provide the expressions for the matching coefficients in the 
limit e — )■ 0, since the 0{e) terms are not relevant in two-loop applications. 

At NNLO the matching coefficients are no longer given by the finite terms of the 
respective form factors alone. We now find 
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The other tensor coefficients are again found to fulfill the four- dimensional constraints 
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which provides a non-trivial cross check of our calculation. 

As a further check of our NNLO results we verified that the matching coefficients 
obey the renormalization group equation, 
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which consists of a universal piece related to the renormalization properties of the SCET 
current with 
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and a second term that contains the anomalous dimension of the QCD current with 
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The one- and two-loop coefficients needed for the check read Fcusp = ACp, t'*-^-* = —5Cf 
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The twofold structure of (jSHD can be used to distinguish the scale /x, that governs the 
renormalization group evolution in SCET, from a second scale u, that is related to the 
non-conservation of the scalar/tensor current in QCD. More explicitly the distinction 
between the scales ii and u can be accounted for by writing 



Cfiu; /i, z/) = Cl{u; /i) + (5C/(m; /i, z/), 



(44) 



where the first term on the right-hand side refers to the above expressions for the matching 
coefficients, C/(u; /i) = Cl{u), while the latter captures the dependence on ln(z//yu), which 
vanishes when the two scales are not distinguished. Expanding the new contribution as 



sa 



E 

k=l 



a 



(4) 



(/i) 



47r 



5C: 



•j,ik) 



we find 



in NLO, and 



5Cp)(«;/i,z/) 



' /x 



(45) 



(46) 



6C: 



Wi), 



u;fx,u) 



7. 



(1)^ 



li Po 



+ 



Ci'^ In^ 






u;fi) 



In- (47) 



in NNLO. (Here /3f^ = 11Ca/3 - 4/3 Tpn/ refers to the QCD beta-function with Uf = 
n/ + 1 fiavours.) Our final results for the matching coefficients with the two scales fi and 
u distinct from each other are provided in electronic form in [23j . 

The matching coefficients with the SCET and QCD scale distinct from each other can 
be used for additional cross-checks. The scalar coefficient is not independent but can be 
related to the vector coefficients by means of the equations of motion, yielding 



Cl{u;fi) 



^2 



>;/^) 



Cy(M;yu) 



rrib 



Cs{u;fi,u) 



(4^ 



where rn}j{v) is the MS renormalized mass in five-flavour QCD. Due to the conservation 
of the vector current the left-hand side of pHj) . which happens to be just the coefficient 
Cf from (IMl) below, is free of v. Hence the QCD scale must also drop out of the 
right-hand side. We checked that our results satisfy ( HHj) . An equivalent formulation 
of fHHl) was given in |9] in terms of a Ward-identity. Also the tensor coefficients at u = 1, 
corresponding to q^ = 0, can be checked against existing results in the literature, since 
they enter the b —> s'j process. From [23] (see also ^25j) one can infer the combinations 



(49) 



- 2 F^{u = 1) + -F|(m = 1) + F|(m = 1) 
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and 

-2C1(m= l;/i,z/)+C3(M = l;/i,z/). (50) 

The latter equation can again be checked for distinct fi and u, and both ( l49l) and ( !50l) 
agree with the formulas in [21]. Note that ( l50l) is just the coefficient Cj. from (15^ at 
M = 1. 

In Figured] we evaluate the matching coefficients for ft = u = rub and as (wife) = 0.22. 
For completeness we show the full set of matching coefficients Cf that we introduced in 
Table [1] We see that the NNLO corrections are in general moderate and add in each 
case constructively to the NLO corrections. In Figure [T] we also show the effect of a 
finite charm quark mass, which is generally rather small, typically modifying the NNLO 
correction by about 10 — 20%. 

4 Exclusive semi-leptonic and radiative B decays 

With the two-loop matching coefficients C^ at hand, we explore several applications to 
B meson decays in this and the following section. For the numerical study we use the 
following input parameters: the 6-quark pole mass rrif, = 4.8 GeV; the renormalization 
scale of the QCD scalar and tensor currents z/ = m^; the hard scale /x = nib. The strong 
coupling constant is obtained from as (rrib) = 0.215 by employing three-loop running 
(A^=^^ = 290.9 MeV), which gives ai^^(1.5GeV) = 0.349. When we add the hard 
spectator-scattering contribution from [H] as required for exclusive processes, we need 
further parameters (such as moments of light-cone distribution amplitudes), for which 
we use the same values as [H] including the hard-coUinear scale /ihc = 1.5 GeV. 

4.1 Heavy-to-light form factor ratios 

The heavy-to-light form factors in the large-recoil regime, where the light meson momen- 
tum is parametrically of order of the heavy-quark mass, take the following factorization 
formula [Hl^ 



Ff^"\E) = aiE)UE)+ / — / dvT,iE;\n 

Jo w Jo 



uj,v)(I)b+{uj)(I)m{v), (51) 



where E denotes the energy of the light meson M, ^a{E) is the single non-perturbative 
form factor (one of two when M is a vector meson), and (px the light-cone distribution 
amplitudes of the B meson and the light meson. The short-distance coefficients Ci and 
the spectator-scattering kernel Tj can be calculated in perturbation theory. The two 
terms in the above equation correspond to the matrix elements of the two terms in the 
operator matching equation ([T]). In particular, the two- loop results from the previous 
section enter the coefficients Ci{E) of the first term. The spectator-scattering kernels Tj 
have been calculated at 0{as) in [7], and at 0{al) in [T3][T^. 

In the following we discuss relations between different QCD form factors Ff^'^^^E) 
that can be deduced from the factorization formula fl3T]) . Adopting the same conventions 
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Figure 1: Matching coefficients C| at the scale ^jl = v = rrn, as a function of u (the 
momentum transfer is given by q^ = {l—u)m^). The dotted horizontal lines show the 
tree level results, the dashed lines the one-loop approximation and the solid lines the 
two-loop approximation with massless charm quarks (orange/light grey) and massive 
charm quarks with nic/mi, = 0.3 (blue/dark grey). 



and notations as [M], we can express the three independent B ^ P form factors as 
f+{E) = df^'\E)^piE) + JdTCf^'\E,T)EpiT,E), 

fo{E) = dfy{E)UE) + JdTC\f'\E,r)Ep{T,E), 



2E 



14 



"^^ ME) = C\f\E)ip{E)+ drCf^'\E,r)Ep{r,E), (52) 

niB + TUp ■'^ J ■'^ 

and the seven independent B ^ V form factors as 

rriB 



my 

■mp + my 
2E 



V{E) = &^'\E)UE) + ldrC\^'\E,T)E^{T,E), 
Ao{E) = &^'\E)i\\{E) + j drCf^'\E,T)^\\{T,E), 
A,{E) = dy^'\E)UE) + 1 drdy'''\E,T)E^{T,E), 



^^^^^ME) - ^^^^^^ME) 
IE nip 



= C}f °) {E) eii {E) + JdT C}f ^) {E, r) S|| (r, E) , 
Ti(E) = d^^'\E)UE)+ [drCS'\E,r)E^{r,E)., 



mp 
2E 



"^UE) = c!,^'\e)UE) + J drC^^'\E,r)E^{r,E) 
T,{E)-n{E) = &;^'\E)i\\{E) + j drC\f\E,r)E\\{r,E). 



(53) 



Here irtp represents the B meson mass, mp and my refer to the pseudoscalar and vector 

hght meson masses, respectively. The coefficient functions Cp and Cp are defined as 

hnear combinations of the matching coefficients of two- ( "AO-type" ) and three-body ( "B- 

type") SCET operators, while '^aij^E) denotes the matrix elements of the three-body 

operators 0\ (si,S2), see ([1]). In terms of the coefficients Cl introduced in previous 

sections, the five independent AO-coefficients are given by 

u 
— ( 
2 



C}^°) = Cl.{u-ii) + '^Cl{u-ii) + Cl{u-ii) 



c 



{AO) 

h 



C^{u- /i) + (l - ^) Cl{u- ft) + C^u; /i) 



4f) = -2C^(«;/i,z/)+(l-^)c|(«;/i,^) + C|(w;/i,^). (54) 

Recall that in D = 4 dimensions one has C|. = C^ = 0. The variable E used in (l52l) and 
fl53|l is related to u through u = 2E/mp. The five independent B-coefficients are given 
in Appendix A2 of [14] • 

From fl52l) and fl53|) . we have the following two identities 

"" ViE) = ^^^^±^ME), T,iE) = ^T,iE) (55) 



mp + mv 2E ' ' ' ' 2E 
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up to power corrections ^E\. In the physical form factor scheme [71[T3], where the SCETi 
form factors ^a{E) are defined in terms of three QCD form factors, 



e/ ^ u, i 



FF 



rriB 



the five remaining form factors read 



V, er - ^^^^7^ ^1 - "^r "^ ^, (56) 



2E 



VflB 



rriB 
2E 



/o — Rq^p + (C'^g — Cj^ Rq] -k Ep, 



rriB 



TUB + rap 



It = i?Ter+(c}^^^-C^i?r)^Hp, 



Ti = R±e/+(c!r':"-cr'R±]^^±, 



,(B1) MBl) 



my 



Ao = i?oer+(c^r-^r^o)^H„, 



rriB 



2E- - ^T^f-ty^f^-^f^iiT)*^\ 



T2-T3 



^(B1)_^(B1) , 



(57) 



In this scheme there are only three non-trivial ratios R and three non-trivial combinations 
of B-coefiicients, defined, respectively, as 
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1 ,"* o(4). 
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V 47r 
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(AO) 



(4) 
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An 
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and 



cr(r.E) 



or^ 



^Mu) 



,{B1) 



47r2 421 



+ Oia^) , (58) 



r,E)-C\:'\r,E)RoiE), 



T,E) 



cSv\r,E) 



C 



a 



/t ^ 
(Bl), 



^,E) 



7+ 
/+ 



r, E) - C, 



(Bl) 
V 

(4) 



(r,E)i?7^(E), 
(r,E)i?x(E). 



(59) 



We denote L^ = ln(/iVm2), L^ = ln(z/Vmg), and /3^^^ = 11/3 Ca - ^/^Tpni. The func- 
tions ji{u) can be found in Appendix [Bl One recognizes the relatively simple structure of 
the ratios Rx in the physical form factor scheme. Compared to the matching coefficients, 
where we encounter up to the fourth power of logarithms, the ratios Rx have logarithmic 
dependences that are at most quadratic, since the universal Sudakov logarithms cancel 
in the ratios. 

As expected in any perturbative QCD calculation, the higher-order correction is nec- 
essary to eliminate scale ambiguities. While the AO-coefficients Cx depend on the 
hard scale /i (which is cancelled by the corresponding dependence of the SCETj form 
factors ^a{.E)), the /i dependence of the ratios Rx {X = 0,T, ±) arises only from the 
scale-dependence of as{fi) and should be reduced after including the higher-order cor- 
rection. In Figure [21 we show the dependence of the three ratios Rx on the scale fi at 
u = 0.85 (corresponding to the light-meson energy E = umB/'i = 2.24 GeV or momen- 
tum transfer g^ = 4.18 GeV^) and fixed renormalization scale u = nih of the QCD tensor 
current. In the absence of radiative and power corrections, all these coefficients equal 
1 (dotted lines). We observe that the scale dependence is reduced at the two-loop order 
for the ratios Ro,t, but not for R±, which receives a large two- loop correction. 

Since the AO-type coefficients C^ and hence the ratios Rx also depend on the 
momentum transfer g^, we show in Figure [31 these coefficients as a function of u (related 
to light-meson energy E = u'mB/2 or momentum transfer g^ = (1 — u)m'^^), with the 
scales fixed at z/ = /i = rrih. As illustrated in Figure [3l the NNLO correction to all the 
five coefficients Cx is quite similar and adds in each case constructively to the NLO 
result; among the three ratios Rx, the two-loop correction to R±, i.e. to the ratio of the 
tensor and vector form factor, Ti/V, is most significant. 

To further investigate these two-loop corrections to the form factor ratios, following 
[Ti] we also take the 5 — )• vr and S — )■ p transitions as examples. Seven ratios among the 
total of ten pion and p meson form factors can be obtained from the two identities f[55l) . 
which do not receive any perturbative corrections, and the five relations that follow from 
flFTl) by dividing through the appropriate ^J^. The g^ dependence of these form factor 
ratios are shown in Figure [H As in f!H] the g^-dependence of the ^^^ in the normalization 
of the spectator-scattering correction is taken from the QCD sum rule calculation. The 
ratios are normalized such that in absence of any radiative and power corrections they 
equal 1 for all g^. Our final results, including both Rx and the spectator-scattering term 
to order a^, are shown as solid dark grey (blue in colour) curves, while the results with 
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Figure 2: Dependence of the ratios Rx (X = 0, T, L) defined in ^5^) on the scale fj,, with 
u = 0.85 (corresponding to the light-meson energy E = wniB/'^ = 2.24 GeF or momentum 
transfer q^ = 4.18 GeV^) and v = rub (the renormalization scale of the QCD tensor current). 
All of them equal 1 in the absence of radiative and power corrections (dotted line). The solid 
and dashed lines denote the NNLO and NLO results, respectively. 



Rx evaluated only at NLO as solid light grey (orange in colour) ones. One can see that 
the radiative correction always enhances the symmetry-breaking effect, and the NNLO 
term is generally quite moderate; the most significant effect from the two-loop correction 
is on the ratio Ti/V (through the ratio Ra_). To see the relative size of the two terms in 
the factorization formula (!5T]) . we also show the result without the spectator-scattering 
term (dashed curves with blue/dark grey and orange/light grey denoting the NNLO and 
NLO results, respectively). Comparing the solid with the dashed curves, one can see 
that the radiative correction from the AO-coefficients Cj^ is always smaller than the 
spectator-scattering contribution. 

To compare our results with the QCD sum rule calculations [2^, the sum rule pre- 
dictions for these form factor ratios are shown as dash-dotted curves in Figure HI One 
notices that, while the sum rule calculation generally satisfies the symmetry relations 
better than predicted on the basis of the heavy-quark limit corrected by radiative and 
spectator-scattering effects, see for instance the lower right panel of Figure HI there are 
also significant differences concerning the sign of the correction, which might be due to 
1/mfe power corrections or ununderstood systematics of the sum rule calculations; further 
detailed discussions could be found in [71[T11|2H]. The new two-loop correction does not 
affect the conclusions on this point. 
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Figure 3: The AO-type coefficients Cj^ and the ratios Rx (X = 0,T, Xj defined in [5^ 
as a function of u (related to light-meson energy E = umB/2 or momentum transfer q^ = 
(1 — u) rn^^), with the scales fixed at v = n = nib- The legend is the same as in Figure\^ 
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Figure 4: Corrections to the B ^ it and B ^>- p form factor ratios as a function of momentum 
transfer q^ . All the ratios equal 1 in the absence of radiative corrections. Solid curves: full 
results with Rx evaluated at NNLO (blue/ dark grey) and NLO (orange/light grey), including the 
spectator- scattering term; Dashed: results without the spectator-scattering contribution; Dash- 
dotted: results from QCD sum rule calculation. The lower right panel shows the two form factor 
ratios that equal 1 at leading power. For comparison, the QCD sum rule results for these two 
ratios are also shown (upper line refers to Ai/V , lower line to T2/T1). 



4.2 Exclusive radiative B decays 

As factorization calculations of exclusive radiative and hadronic B decays involving only 
light mesons make use of the heavy-to-light form factors at maximal recoil, it is of interest 
to investigate the short-distance corrections at m = 1, i.e. E = mB/^ or q^ = 0. In this 
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subsection we shall consider the following two ratios 



7^l(^) 



7^2(^) 



TTiB + mpUiE) ^^ ' Jo ^+ ^ ' ^ /+(E) 



niB + mvTi{E) 



R±{E) 



niB + niv l'\,_^{Bi)^^ j^^^±ir,E) 



[ '-^nv(r,E)^^ 



TUB y{E) ■ ■ niB 

defined in the physical form factor scheme. 

At u = 1 and assuming the asymptotic form for the light-meson distribution amplitude 
4>m{v) = Qvv, the analytic expressions for these two ratios simplify considerably, even 
at NNLO. As the spectator-scattering contribution is already given by Eq. (124) in [13], 
here we give only the expressions for the ratios Rt,± at m = 1 (as a consequence of the 
equations of motion, we have Ro{u = 1) = 1), 
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(61) 



with L^ = ln(yU^/m^), Ly = ln(z/^/m^), and ni = A has been used. Using the three-loop 
running coupling and specifying to the pion (T^i) and p meson (7^2)5 numerically we 
obtain (setting u = fi = mi,) 



TZiiE, 



1 + 



0.046 (NLO) + 0.015 (NNLO) {Rt 



-0.160 |l + 0.524 (NLO spec.) - 0.002 {5l^] 



0.817, 



T^2{Emax) 



1- 



0.023 (NLO) + 0.030 (NNLO) (i?. 



+0.084 <^ 1 + 0.406 (NLO spec.) + 0.032 (S 



logy 
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(62) 



In these expressions we separated the symmetry-conserving (first number, normalized 
to 1), AO- and B-type corrections (denoted by Rt,± and the remaining terms, respec- 
tively). The parameter Sf denotes the small effect from renormalization-group summa- 
tion and has the same meaning as in Eq. (124) of ^14j. We observe that the AO-type and 
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spectator-scattering corrections always have opposite sign, but the latter are larger and 
determine the sign of the deviation from the symmetry limit. We also notice that the 
two-loop correction to R± is more significant than to Rt- The small numerical difference 
of spectator-scattering contribution relative to Eq. (124) in [H] is due to the fact that now 
the 3-loop running coupling is used. For comparison the QCD sum rule calculation [27] 
gives 7^1 = 0.955 and 7^2 = 0.947. For the tensor-to- vector ratio 7^.2, one notices that the 
sign of the symmetry-breaking correction between these two methods is opposite. Since 
the form factor ratio Ti/V is important for radiative and electroweak penguin decays (see 
the discussion in Section 5.2 of [M]), the discrepancy between the SCET and QCD sum 
rules results for 7^2 suggests that a dedicated analysis of symmetry breaking corrections 
to form factors (rather than the form factors themselves) with the QCD sum rule method 
should be performed. 

5 Semi-inclusive B — > Xsi^i~ decays 

Rare inclusive 5-meson decays induced by the quark level transition b — )■ si~^i~ are highly 
sensitive to new physics. Due to the presence of two extra operators {ii)v,Ai^^)v-A in 
the effective Hamiltonian and the availability of additional kinematical observables, such 
as the dilepton invariant mass (g^) spectrum and the forward-backward asymmetry, the 
b — )■ si~^i~ decay provides complementary information relative to the radiative b ^ s'j 
process. 

The exclusive decay process B — > K*i~^i~ has been studied in great detail, both 
with respect to its QCD dynamics [29j and to the sensitivity of various observables 
to new physics [30], because it can be measured relatively easily at hadron colliders. 
Also on the inclusive decay process B — )■ Xsi~^i~ dedicated work exists on higher order 
radiative corrections (see [21] for recent reviews), power corrections [521E3], and on the 
identification of additional kinematic observables [34j. 

The low dilepton invariant mass region, 1 GeV^ < Q'^ < 6 GeV^ is particularly in- 
teresting, since it benefits from smaller theoretical uncertainties and a higher rate. At 
somewhat higher q"^ the spectrum is dominated by charmonium resonances (which also 
determine the integrated decay rate, see the discussion in [53]). On the other hand, for 
q^ < 1 GeV^, the branching ratio is determined largely by the contribution from almost 
real intermediate photons, and hence contains essentially the same information as the 
6 — 7- S7 transition. 

In the following we discuss semi-inclusive B -^ Xsi^i^ decay, where the hadronic 
final state Xg is constrained to have small invariant mass rrix and g^ is in the range 
from 1 GeV^ to 6 GeV^. In this kinematic region (the so-called "shape function region"), 
the outgoing hadronic state is jet-like and the relevant degrees of freedom are hard- 
coUinear and soft modes. The semi-inclusive decay rates can be calculated by matching 
the effective weak interaction Hamiltonian to soft-coUinear effective theory. At the leading 
order in the ^.qcul'^h expansion, the decay rates can be factorized into pro cess- dependent 
hard functions h^^\ related to physics at the hard scale /i ~ rrife and above, a universal jet 
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function J, related to physics at the intermediate hard-colhnear scale /ihc ~ \nniAQCD, 
as well as a universal non-perturbative shape function S, describing the internal soft 
dynamics of the B meson, with the following schematic form 



drM = /iM X J0S, (63) 

a result already applied extensively to inclusive B — )• X„£z/ and B — )■ Xs'j decays in 
the shape-function region. The two-loop matching coefficients of the tensor currents 
calculated in the present paper provide further input to reaching NNLO (a^) accuracy in 
/if°l and the entire differential decay rate dr'*''. Compared to exclusive decays mediated 
by the b — )■ si~^i~ transition [29] the semi-inclusive case has the advantage that the 
theoretically less certain spectator-scattering contributions to the currents that enter the 
exclusive form factors are power-suppressed and can be dropped. 

In the following we will be mainly interested in the forward-backward asymmetry of 
the differential rate integrated up to an invariant mass m^x^ in the final state. We briefiy 
review the theoretical description of this quantity, adopting the same conventions and 
notation as |37], to which we also refer for further details. The short- distance coefficients 
/it"! at the hard matching scale fx are composed of products of two factors, since the 
hadronic part of the effective weak interaction Hamiltonian is first matched to two QCD 
(rather than SCET) currents. 



4 = s rPLb , J!; = '^rs iq^a^^PRh 



(64) 



u=mf, 



with coefficients Cf^^^{q^,jj) and Pl,r = (1 =F 75)/2. Moreover, m^ in J^ refers to the 
bottom quark pole mass. The QCD currents are then related to the corresponding SCET 
currents. 



4 = E ci{u,^,)[lWr.,]T^,^,K. 



j=l,2,3 
^ i=l,2 

These equations represent the momentum space versions of ([1]). The variable u is related 
to the kinematics of the process by u = p~ /mb, where 

p~ = n+p = rrib -r, (66) 

niB -Px 

and Px = n_px ^ ttib is the small light-cone component of the hadronic final state's 
momentum. The basis of Dirac structures is chosen as 
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As noted in ^7\, the choice of qf^ instead of nt for Fgg is convenient here as it makes 
exphcit the constraint from lepton current conservation, which imphes that for massless 
leptons Cg does not contribute, while for Ti^- there are only two independent coefficients. 
Transforming the basis ( 167|) to our operator basis listed in Table [H the matching coeffi- 
cients c^ and c[ are given, respectively, as 



cliu^ 


./^) = 


CUu; 


/^), 








cl{u. 


./^) = 


cliu- 


^l) + 


-CUu 
u 


;/^), 




cliu. 


,/^) = 


2 

umi 


-cl{u■^l), 






cliu. 


,/^) = 


:-2C^(M;/i, 


u = mt 


,) + C^iu: 


/i,Z/ 


cl{u. 


>/^) = 


2 


-cUv 


t;/i, i/ = 


rrib). 





rrib) , 

(68) 

urrih 

The two-loop matching coefficients c^ for the vector current have become available in the 
context of inclusive semi-leptonic B decays P4T2]. The results of this paper allow us to 
compute also the matching coefficients cj at NNLO. As a consequence the factor in h^'^^ 
related to the QCD current matching is now complete at NNLO, while the other factor 
related to C™^^(g^,/i) is known at the next-to- next-to-leading logarithmic (NNLL) order, 
since the three- loop 0(0^) matrix elements of the current-current operators (giving rise 
to charm- loop diagrams) are not available. 

In Figure [5] we show these matching coefficients as a function of u in the one- (dashed) 
and two-loop (solid) approximation, evaluated at /i = mft = 4.8 GeV (blue/dark grey 
curves) and aX ^ = 1.5 GeV (orange/light grey curves), respectively. The difference 
between these two different choices of the IR factorization scale n is compensated by 
the corresponding scale dependence of the convolution J ^ S such that the differential 
rate f l^ is //-independent. Note that, while we show the entire range of u, Eq. f pB]) 
implies that the relevant values of u for b — )■ si~^i~ in the g^ region of interest are above 
u ~ 0.75. In the lower right panel of Figure |5l we also show the ratio cl/cf, which equals 
the quantity R± defined earlier in fl58l) at u = rrih, and plays an important role for the 
forward-backward asymmetry as discussed below. Note that R± is /i- independent, except 
for the truncation of the perturbative series. In evaluating this ratio to a given order in 
as, we expand the denominator and truncate the expanded expression. 

Comparing the dashed (one-loop approximation) and solid (two-loop approximation) 
curves of the same colour in Figure [5l we observe that the two-loop corrections are 
generally moderate in the large u (low g^) region, whereas the large correction in the 
region of small u is due to the fact that increasing powers of large logarithms take over 
in this region. However, the correction is amplified in the ratio R±, where the two- loop 
correction exceeds the one-loop term. This leads to a considerable residual /i-dependence 
(difference of blue/dark grey and orange/light grey curves) as can also be seen in Figure [2j 
Since the infrared physics drops out from the ratio c[/cf the natural scale is of order of 
the hard scale m^. 

24 




0.0 
-0.1 
-0.2 



-0.3 -, 



-0.4 L^ 
0.0 




0.2 0.4 0.6 
u 



0.8 1.0 



0.15 
_. 0.10 

3. 

^ 0.05 



0.00 




0.0 0.2 0.4 0.6 
u 



0.8 1.0 




Figure 5: The matching coefficients c^{u,^) and cj{u,^) as a function of u (related to the 
dilepton invariant mass g^ = (1 — u)m1) in the one-loop (dashed) and two-loop (solid) approxi- 
mation. The blue/dark grey curves refer to ^ = nif, = 4.8 GeV, and the orange/light grey ones 
to ^1 = 1.5 GeV. 



The differential decay rate (l63l) can be written as 



d^T 



dq^dpxd cos 6 8 

+ 2 cos eHA{q\p\) 



(1 + cos" d)HT{q\p\) + 2 (1 - cos^ d)HL{q\p\) 



(69) 



where for B decay, 9 denotes the angle between the positively charged lepton and the B 
meson in the centre-of-mass frame of the £+£~ pair. For fixed pj^, the forward-backward 
asymmetry in 9 therefore vanishes for a particular q^ at which HA{qQ, Px) = 0- Integrating 
over the invariant mass of the hadronic final state up to the cut m3^*, the asymmetry 
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zero occurs at 



■Pi 



dp\HA{ql,p\) 



P+-* 



const X / - dp\h^'X^lp\)^^^^-^^^l I dujp-J{p~u)S{p%-uj), (70) 



9o^ 



where 



q+ = mB-Px, q- 



q^/q^ 



+cut 

Px 



2m 



B 



ml + {m^y - g' - 



^(m| + (mT) 



cut\2 



^2^2 



4mn(m3^*)2 



•'B\ 



and 



/iA(g',Pi)=2Cioc?(M)Re 



cr'(g')c?(M) + — ^^'^'^^^"^^ 



Cr^(g2)cl(«) 



(71) 



(72) 



We now observe that ^V(q'o)Px) depends on p^ only through the definition of u in 
( l66l) and the kinematic factor 2mh/q-. For typical r7i3f* of 2 GeV this dependence is very 
weak, since then pj^ ~ 1 GeV <^ m^. Thus, Px appears only as a small correction to 
ruB — Pxj cind in the definition of m in a term that is additionally suppressed by q^ /mB 
relative to m;,, see fl66l) . This results in a very small variation of u of about 0.02 over 
the entire pj^ integration region. We may therefore pull the slowly varying function 
^AWO'^x) ^^ front of the pj^ integration in (I70p thereby replacing pj^ in the argument 
by an average value which we assume to be (pj^) = pJ^^^*/2. The remaining integral over 
the jet and soft function is different from zero, thus the forward-backward asymmetry 
zero is determined by hj^{ql, {Px)) = 0. Using (172!) this is equivalent to the condition 



ll 



2mb{mB - {Px)) 



Re[Cr'(go')]c?(%) 



(73) 



with uq = 1 — qQ/{m,b{m,B — (Px)))- This result leads to the important conclusion that 
the QCD dynamics that determines the location of the asymmetry zero is to a very good 
approximation independent of the long-distance physics below the scale mi, contained in 
the jet function and the non-perturbative shape function. It also depends only very weakly 
on the value of the invariant mass cut through the dependence of (pj^) on m^^*. The bulk 
dependence of q^ on the invariant mass cut m™* enters through the kinematical factor 
ruB — (Px) on the left-hand side of fl73|) . 

We are now in the position to quantify the impact of the two-loop calculation of 
R±{uo, V = nih) = cl{uo)/cl{uo) on q^. In ^7\ the asymmetry zero has been determined 
by keeping the full NNLL expression for Re [Ci^''''\q^)]/Re [C'^''\q^)] but setting R^ = 1. 
In this approximation, and excluding l/r/ib-suppressed shape function effects for the 
moment, the zero is found to be 



^h 



(3.62 . . . 3.69) GeV^ for 



m 



X 



(2.0... 1.8) GeV. 



(74) 
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as{Mz) = 0.1180 


A2^i(m|.-m|)~ 0.12 GeV' 


sin^ Ow = 0.23122 


^polc _ ^^^_^ Q^y 


Mw = 80.426 GeV 


mP°i^ = (1.5 ±0.1) GeV 


Mz = 91.1876 GeV 


mPS(2 GeV) = (4.6 ± 0.1) GeV 



Table 2: Numerical inputs that we use in the phenomenological analysis of the forward-backward 
asymmetry zero. 



As indicated the lowest value corresponds to m™* = 2.0 GeV and the highest one to 
772™* = 1.8 GeV. Our value is somewhat larger than what can be extracted from Figure 4 
of ^37j, because we expand the factor mb(/i)/m^°'^ that accompanies Cj in as- Moreover, 
the variation of the zero when changing m™* from 1.8 GeV to 2.0 GeV is about twice 
as large compared to what can be read off from Fig. 4 of [37] , which is likely due to our 
approximation of pulling the slowly varying function ^^4 (O'o'Px) "^^^ ^^ ^^^ integral in 
(!70|) . However, our approximation is still justified since even the increased sensitivity of 
the zero on m™* is only ±0.03 GeV' and hence below 1%. Taking into account -R_l at 
the NLO, we find for the position of the zero 



<i\ 



(3.55 . . . 3.61) GeV' for 



m 



cut 
X 



(2.0... 1.8) GeV. 



(75) 



I R.^ NLO 



The impact of the NLO correction to i?_L is to shift the zero by —2.2%. As we already 
stated before, and as can also be seen from Figures [2] and 0, the size of the NNLO 
correction to R± is significant. It amounts to a shift of the NLO zero in fl75l) by another 
—3% and hence is larger than the NLO shift. The total shift induced by R±_ through 
NNLO therefore amounts to —5%. 

Before proceeding to our final result we briefly comment on the role of power correc- 
tions. The authors of [27] performed a thorough study of l/mf,-suppressed shape function 
effects which result in a shift of the zero of —0.05 GeV' to —0.1 GeV'. This shift is more 
strongly dependent on the invariant mass cut and the theoretical error increases when 
m™* is chosen smaller. In the following we take the larger value as an estimate for the 
shift and also for the associated uncertainty. However, the study of power corrections 
in [37] does not cover all such corrections and applies a rather crude treatment to those 
arising from soft gluon attachments to the charm-loop diagrams by absorbing the 1/m' 
non-perturbative power corrections into the C™^\ which is justified only in the absence 
of invariant mass cuts. In the semi-inclusive region, the matrix element of (29) in [33] 
cannot, due to the presence of a soft gluon, be expressed in terms of a short-distance coef- 
ficient times a local matrix element, since the soft gluon attached to the charm loop affects 
the invariant mass of an energetic hadronic final state by a relevant amount -^mftAQCD, 
which must be accounted for by a subleading shape function. By treating this correction 
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as in the inclusive case, the authors of [2Z] imphcitly assumed that this shape function 
somehow factorizes into the local heavy-quark effective theory matrix element A2 and 
the leading-power shape function. It is not clear to us how this simplification can be 
justified and it is likely not even parametrically correct. Nevertheless, in the absence of 
better information we follow the treatment of [37] and include the l/rn^ power correc- 
tions into the Q"'^'. This results in a shift of the asymmetry zero by +0.07 GeV^, which is 
included in ( !71|) . (175|) . and below in (176|l . To be conservative we assign another 0.1 GeV^ 
uncertainty to this estimate and add it in quadrature with the other power correction 
uncertainty. 

We are now in the position to present our final NNLO result based on the numerical 
input parameters and their respective intervals as specified in Table HI We then find 

q^ = [(3.34... 3.40) tHa ^ ± 0.08 „,+H4n^.±0. 14 sF±0. 14 ^p+^]GeV2 

= [(3.34 . . . 3.40)+°:22] GeV^ for mf = (2.0 . . . 1.8) GeV . (76) 

The error estimate is computed as follows: The range of scale variation is taken to be 
2.3 GeV < /i < 9.2 GeV, and we vary the scale in the Cf"^^ and in R^ independently to 
account conservatively for the absence of the 0{al) correction to the C™'^^ The input 
quark mass is the bottom mass in the potential-subtracted (PS) scheme |38], see Table |2j 
The pole mass and MS mass used in intermediate expressions are computed using the 
one-loop conversion factors resulting in m\°'^ = 4.78 GeV and m{rn^'^) = 4.36 GeV, 
respectively, when m^^(2GeV) = 4.6 GeV. The dependence on the charm quark mass 
enters through the matrix elements of the current- current operators. The error labelled 
"SF" is connected with the subleading shape function effects as discussed above. Finally 
we have added an uncertainty estimate for the approximation made by pulling out the 
slowly varying function /1I4 (^oiPx) out of the pj^ integral in (I7U]) . We estimate this error 
rather generously by varying (pj^) from p3f*/4 to 3p3f*/4. The total error is obtained by 
adding all these uncertainties in quadrature. 

We note that the value of the asymmetry zero in semi-inclusive b — )■ s£"*"£~ decay 
is significantly smaller than for the exclusive case [29|, where spectator scattering is 
responsible for a positive shift as is the fact that in this case {p~^) = in fl75]) . On the 
other hand the semi-inclusive zero is in the same region as in the inclusive case [39lj . 
where virtual effects together with hard gluon bremsstrahlung encoded in functions uno 
and cugio [^ also induce a negative shift on the zero. 

6 Conclusion 

In this paper we completed the two-loop matching calculation for heavy-to-light currents 
from QCD onto SCET for the complete set of Dirac structures. These matching co- 
efficients enter several phenomenological applications, of which we have discussed their 
effects on heavy-to-light form factor ratios, exclusive radiative and semi-leptonic decays, 
as well as the inclusive decay B — )■ Xsi~^i~ in the shape-function region. The two-loop 



28 



corrections are generally relatively small, in the few percent range. However, one ratio, 
R± = c[{u,fi)/cl{u,fi), which is also the most important for phenomenology, since it en- 
ters the comparison of radiative and semi-leptonic decays as well as the forward-backward 
asymmetry in exclusive and semi-inclusive b — )■ si~^i~ transition, exhibits a two-loop cor- 
rection that is larger than the one-loop term. The two-loop term alone shifts the location 
of the asjTumetry zero by about —0.1 GeV^, comparable to the effect of l/rrif, suppressed 
shape functions estimated in (3^. We showed that the location of the asymmetry zero in 
semi-inclusive B -^ Xsi~^£~ with an invariant mass cut is to a very good approximation 
independent of the long-distance physics below the scale rrib contained in the jet function 
and the non-perturbative shape function, and obtain Qq = (S.SA^q'H) GeV^ for an invari- 
ant mass cut m™* = 2.0 GeV as our best estimate for the asymmetry zero. Moreover, 
we confirm the discrepancy between QCD sum rule and SCET results for the form fac- 
tor ratio Ti/V in the low g^ region discussed in |T1] and suggest that a dedicated QCD 
sum rules analysis of deviations from the symmetry limit (rather than the form factors 
themselves) should be done to clarify the situation. 
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A NLO coefficient functions 

In Section 13.11 we introduced the following set of one-loop coefficient functions, 

5 
9oiu) = -- + 21n(M), 

71-2 2 

gi{u) = -— + - \n{u) - 2 \n\u) - 2Li2(u), 
12 u 

/ N ^^ 1^ 12(l+M)+7r2M 2/, 2/ N T- / nA 4,3,, 

92{u) = ^ + 3C3 + ^ ^ ln(M) - - (^ \n\u) + U2{u)) + - \n\u) 

+ 4 ln(u)Li2(u) - 2Li3(u) + 4Si,2(m), 

^'^""^ = "160 " 6^=^ + ^ 12^ ^ ^''^''^ ~ 3 ^""^""^ ~ ^ ln'(i^)Li2(u) 
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5-4 (m) = gi{u) - 6 - — ln{u) 
u 



95{u) = g2{u) - 10 - — - — 31n(M) -2\n\u) - 2U2{u)), 

9,{u) = g,{u) - 18 - ^ + 1(3 - - f^^^ Hn) - 3 In^lix) - ?,U,{u) + ^ ln=^(«) 
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u\ 3 

fl'l2(M) =5'7(m) + 2, 

2m 
^i3(m)=^8H + 6 + — ln(M). (77) 

B NNLO coefficient functions 

The finite parts of tlie two-loop form factors involve the following coefficient functions, 
2(7-2m + 3m^ 4(11 + 2iz + 3iz2) qwuw-^ 
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2m 
2m(1 + 3m) 



2ho(u) 



u-^ 



8S2,2(M)+2Li4(M)-Li2(M 



„-^\2 



20 



2(3 + 21m - 24m2 + 2m3) 4(1 + 7m + 4m2) , , ^ , , 

—^ Li3(m) ^3 ln(M)Ll2(M) 

2(1 + 69m - 48m2 - 24m3 



3m3 



3m2 

o / N / . 2(l-llM-2M2)7r2\ , , , 

■Si,2(«) + ( 4 + ^ — '— ] ln(M) 



32 



2(9 - 13m - 18^2) 2u(1 + 3u)tt^\^ , , 7+15m, o. ^ 



3m^ u^ J 3u 



— ( Li3(-n) - ln(M)Li2(-M) ln(l + u) 

{11 + 17u + 38u^)tt^ 4(3 -u)^ ^i ..n 5435 

^ ^ + ^ ^Cs - 4t^ ln(2) 



9m2 m ^•^ ^ ^ 108 ' 

. / N 1, / N 1, . A 87r2 181 

Moreover, for the ratios Rx in fISS]) we need the following auxihary functions, 

, , 4(m - 2) (u^ + 2m - 2) , , 16u , , 8(2u - 7)u , , 4u , , 

2{u + 3){u^-l) 16(2M-3)MLi4(M) A{u'^ - 3Qu + 2h)U2{u) 



+ ^ S5(m)- 

■U'^ U"^ w^ 

4 (4^3 - 8^2 + 5m + 3) (Li3(M) - C(3)) 4 (Tm^ + 85^^ - niu + 3) Li3(M) 

87r2(6M-17)MLi2(n) 6nLi2(M) 4 (5^3 + 63^2 - 83u + 3) Li2(n) ln(u) 

3'U'^ u u^ 

27r^(20M - 33)m 8 (m^ + 51m - 62) C(3) 271^ (9^2 _ 73^ + 59) 

45m^ v? 3v? 

327r2(5M - 6) ln(M) 2 {3u^ + 41^^ - 55^ + 3) \xi^(u) \n{u) 

4(2M-l)2ln2(M) (9M-4)ln(M) 2(11m - 25) In^H 

. , , 2(m - 2) (m^ + 2m - 2) , , 8m , , 2(4m - 7)m , , 

J2(m) = ^-^z Si(m) - -^ S2(m) + —- S3(m) 

3u'^u 3u-^ 3u-^ 

{u + 3){u^-l) , , 27r2(12M-29)uLi2(M) 4(4m - 15)MLi4(u) 
sgfw) + 



M^ 3u^ V? 

2 (25'u2 - 99m + 51) Li2(M) 2 (3^3 + 29^^ - 37m + 3) (Li3(M) - C(3)) 

3u^ u^ 

2 (11m3 - 63^2 + 57m + 3) Li3(M) 2 {Su^ - 48m2 + 43m + 3) Li2(M) ln(M) 

M^ M^ 

22Li2(M) 47r^(5M-18)M 2 (7m2 - 83m + 86) C(3) vr^ (8m2 - 69m + 67) 
3m 45m^ m^ 3m2 

47r2 (2m2 - 17m + 18) ln(M) (5m3 - 33m2 + 29m + 3) ln2(M) ln(M) 131n2(M) 
3m2 M'^ 3m 

47r2ln(M) 2031n(M) (8m - 51) In^(M) 2571n(M) 269 
3^^ ^ 9^^ 3m 9 ^~9~' 
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■ / N 26 , , 8"" ri 2/ \ T • /-\i 76 

Jsiu) = -— gg{u) + — [In [u) + Li2(m)J - — , 

327r^{u + 2)u 32M(Li3(M) - C(3)) 8(m - 2) (m^ + 2m - 2) Lialw) 

52u\n{u) 104 801n(M) 236 
9^ 3u 3u 9 ' 

4(m - 2) (m2 _ 2m + 2) ,, 16m , , 8(2^2 - 12m + 11) ,, 

4 , , 2(m + 1) (m2 + 2m + 7) , , 16(2M2-8M + 7)Li4(M) 

+ ;— S4{U) -— S5(m) + 



+ 



3m 3m'^ m^ 

2 (51m2 - 328m + 250) U2{u) _ 4 (m^ + 48m2 - 69m - 7) (Li3(M) - C(3)) 
3m2 3m3 

4(3M3 + 315M2-519M + 7)Li3(M) 87r2(6M2 - 32m + 29)Li2(M) 2Li2(M) 



3m^ 3m'^ u 



4 (2m3 - 237m2 + 387m - 7) LIsIm) ln(M) 27r4(20M2 - 83m + 73) 

3m'^ 45m^ 

8(4m2 - 183m + 306)C(3) 7r2(llM2 - 206m + 218) 87r2(M2 - 72m + 120) ln(M) 



3m2 3m2 9m2 

2(7M3-159M2 + 255M-7)ln2(M)ln(M) 41n2(M) ln(M) ,,, 

+ ^ TT-. — — -^ + r^^ + -^ - 13 ln(M) 

3m'^ u u 

2(27M-125)ln2(M) . , , 563 

- ^ ^^ - Sn^ n 2 + , 

3m ^ ^ 24 ' 

.,, 2(m - 2) (m2 - 2m + 2) ,, 8m ,, 2(7m2 - 17m + 11) ,, 
M-) = ^^ -i(-) - ^ ^^(-) ^s -3(.) 

(m + 1)(m2 + 2m + 7) , , 4(7M2-25M+19)Li4(M) 



H TT^ S5[U) + 



3u^ u^ 

2 (35m2 - 133m + 73) U2{u) 2 (7m3 - 99m2 + 129m - 7) (Li3(M) - C(3)) 

3u^ 3u^ 

2(45M3-273M2 + 273M + 7)Li3(M) 27r2(2lM2 - 59m + 41)Li2(M) 22Li2(M) 



+ 



+ 



3m3 3m3 3m 

2 (28m3 - 204m2 + 207m + 7) Li2(M) ln(M) 7r4(35M2 - 122m + 92) 
3m3 45m'^ 



2(40m2-369m + 378)C(3) 7r2(68M2 - 279m + 267) 47r2(8M2 - 75m + 78)ln(M) 



+ 



3m2 9m2 9m2 

(llM^- 135M2 + 14lM + 7)ln2(M)ln(M) 13 ^^(m) 2691n(M) An'^\n{u) 



3m3 3m 9m 3m 



215 In M 40M-73 In^ M , o, / x 4421 

+ ^ - ^ '- ^ + A-K^ In 2 

9 3m ^ ' 216 



34 



. , , 38 , , Su ... . ^ . , ,, 47r2 205 
JrW = — ggiu) - — [In [u) + Li2(m) 



9 "^ ' 3m L ' ' -V 'J 9 54 ' 

87r2(M2 + 8M-16) 32n(Li3(u)-C(3)) 8(m - 2) (n^ - lOn + 10) Li2(n) 

761n(M) 232 4(19m - 156) ln(M) 1429 
^ 9m ^ 9m ^^F' 

57r2(5M + 4) 16M(Li3(M) -C(3)) 2(m - 2)Li2(M) u^\n^{u) 



J9W = — —^ + _2 

563 



6m m m m 

- 12 ln(M) - 6C(3) + An^ ln(2) + 



16 ' 

, , 47r2(M + l) 8M(Li3(M) - C(3)) , , , ,, , o, , , 5141 , , 

^''^""^ = 3^ + J2 + 4 ln(«) + 3C(3) - 27r2 ln(2) - ^^ , (79) 



with 



si(m) = UHiiu) + vr^ ln(2 - m) , 

S2(m) = 12H2(m) - TT^Lisl-M) , 

ITtt^ 



S3(M)=3Li^(M)-24S2,2(^) ^^ , 

60 

s^{u) = QU^iu) - 3Li2(M) ln(M) + 3Li3(M) - 27t^ ln{u) - 6C(3) , 

S5(m) = -2Li3(-M) + 2Li2(-M) ln(M) + ln(M + 1) ln^{u) + n^ ln{u + 1) . (80) 
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